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ABSTRACT 

We derive lens distortion and magnification profiles of four well known clusters observed with Subaru. Each 
cluster is very well fitted by the general form predicted for Cold Dark Matter (CDM) dominated halos, with 
good consistency found between the independent distortion and magnification measurements. The inferred 
level of mass concentration is surprisingly high, 8 < Cyh < 15 ((cvh} = 10.39 ±0.91), compared to the rela- 
tively shallow profiles predicted by the ACDM model, = 5.06 ±1.10 (for (Mvir) = 1.25 x 10'^Mq//i). This 
represents a 4a discrepancy, and includes the relatively modest effects of projection bias and profile evolution 
derived from N-body simulations, which oppose each other with little residual effect. In the context of CDM 
based cosmologies, this discrepancy impUes clusters collapse earlier (z > 1) than predicted (z < 0.5), when the 
Universe was correspondingly denser 

Subject headings: cosmology: observations - gravitational lensing - large-scale structure of universe - galax- 
ies: clusters: individual (A1689, A1703, A370, RXJ1347-11) 



1. INTRODUCTION 

Dark matter (DM) is understood to comprise ~ 85% of the 
material Universe (Fukugita & Peebles 2004), in the form of 
massive halos around galaxies and clusters of galaxies. The 
density profile of a dark halo depends on the unknown nature 
of the DM and the way structure develops over cosmic time. 
Galaxy halos may be modified significantly when gas associ- 
ated with the disk of a galaxy cools and condenses during its 
formation. In contrast, clusters are so massive that the virial 
temperature of the gas (typically, 3-15 keV) is too high for 
efficient cooling and hence the cluster potential reflects the 
dominant DM. 

Models for the development of structure now rest on ac- 
curate measurements of the power spectrum of mass fluctua- 
tions, the cosmological density of DM, and the contribution 
of a cosmological constant, A (Spergel et al. 2007; Komatsu 
et al. 2008). This framework is the standard "ACDM" cosmo- 
logical model, with the added simple assumptions that DM is 
colUsionless, reacts only gravitationally, was never relativis- 
tic ('Cold'), and with initially Gaussian-distributed density 
perturbations. In this context, detailed A-body simulations 
have established a clear prediction that CDM-dominated clus- 
ter halos should have relatively shallow, low-concentration 
mass profiles, where the logarithmic gradient flattens contin- 
uously toward the center with a central slope tending towards 
r"', interior to a characteristic radius, < 100-200kpc/i"' 
(Navarro, Frenk, & White 1997, NEW, hereafter; Bullock et 
al. 2001; Maccio et al 2006; Hennawi et al. 2007; Neto et al. 
2007; Duffy et al. 2008). This "NEW" profile is characterized 
by the (total) mass, Mvii , within the virial radius, rvii , and by 
the concentration parameter, Cyir = rvu/rs- 
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Interaction between clusters indicates DM is collisionless, 
in particular the "bullet cluster"(Markevitch et al. 2002) - 
where shocked gas lies between two substantial galaxy clus- 
ters, implying these clusters passed through each other re- 
cently (Markevitch et al. 2002). Here the weak lensing signal 
follows the bimodal distribution of member galaxies reflect- 
ing substantial DM associated with each cluster component 
(Clowe et al. 2006) and that this DM is relatively collision- 
less like the galaxies. These observations disfavour the class 
of alternative gravity theories for which the lensing signal is 
expected to trace the dominant baryonic contribution of the 
gas (Clowe et al. 2006). Other cases of interaction show that 
in general displacement of the gas relative to the DM is typi- 
cally related to interaction (Jee et al. 2006; Okabe & Umetsu 
2008). 

Eor many clusters no obvious evidence of recent interac- 
tion is seen as the gas and member galaxies follow a sym- 
metric, structureless distribution. Measurements of the mass 
profiles of these relaxed clusters may help in understanding 
the nature of DM, preferably relying on gravitational lensing 
signals where model-dependent assumptions are not required. 
For A 1689, over 100 multiply-lensed images have been used 
to derive the inner mass profile (Broadhurst et al. 2005a), with 
the outer profile determined from weak lensing (Broadhurst et 
al. 2005b). Together, the full profile has the predicted NEW 
form, but with a surprisingly high concentration when com- 
pared to the shallow profiles of the standard ACDM model 
(Broadhurst et al. 2005b; Umetsu & Broadhurst 2008). Fur- 
thermore, good consistency is also found between the lensing 
based mass profile and the X-ray and dynamical structure of 
this cluster in model independent analyses (Lemze et al. 2007; 
Lemze et al. 2008, in preparation). A similar discrepancy is 
also indicated by lensing observations of other clusters e.g., 
MS2137-23 (Gavazzi et al. 2003) and CL0024H-16 (Kneib et 
al. 2003). 

Here we examine a sample of relaxed, high-mass clusters, 
to test the distinctive prediction of ACDM, and to settle em- 
pirically whether the mass profile of A1689 is unusual. In 
section § 1 we describe the data reduction in §3 we present the 
weak lensing analysis and the magnification profiles derived 
from background number counts and make comparison with 
the prediction of the ACDM; in §4 we discuss our results and 
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conclusions. 

2. DATA REDUCTION 

We analyse deep images of A1689, A1703, A370, 
RXJ1347-11 taken by the wide-field camera Suprime-Cam 
(34/ X 27/, Miyazaki et al. 2002) in Subaru (8.3m) which are 
observed deeply in several optical passbands, listed in Table 1 . 
These clusters are of interest due to the exceptional quaUty of 
the data, with exposures in the range 2000s- 10000s per pass 
band, with seeing ranging from 0".5 to 0".75. We use ei- 
ther the R or I bands for our weak lensing measurements (de- 
scribed below §3) for which the instrumental response, sky 
background and seeing conspire to provide the best images. 
All the available bands are used to define colours with which 
we separate cluster members from the background in order 
to minimise dilution of the weak lensing signal by unlensed 
objects. 

The standard pipeline reduction software for Suprime-Cam 
(Yagi et al. 2002) is applied for flat-fielding, instrumental 
distortion correction, differential refraction, PSF matching, 
sky subtraction and stacking. Photometric catalogs are con- 
structed from stacked and matched images using Sextractor 
(Bertin & Arnaut 1996). We select red galaxies with colours 
redder than the colour-magnitude sequence of cluster E/SO 
galaxies. The sequence forms a well defined line due to the 
richness and relatively low redshifts of our clusters. These 
red galaxies are expected to lie in the background by virtue 
of k-corrections which are greater than for the red cluster se- 
quence galaxies and convincingly demonstrated spectroscop- 
ically by Rines & GeUer (2008). We also include very blue 
galaxies falling far from the cluster sequence,> l.'"5-2."'5, 
depending on the passbands available for each cluster (listed 
in Table 1), to minimise cluster contamination (see Medzinski 
et al. 2007). Typically the proportion of blue galaxies used is 
around 50% of the red background. In addition, we adopt a 
conservative magnitude limit of m < 25.5-26.0, depending 
on the depth of the data for each cluster, to avoid incomplete- 
ness. Imposing these strict limits is a necessary precaution 
against the diluting effect that unlensed cluster and foreground 
galaxies will otherwise produce, leading to spuriously shallow 
central mass profiles, as shown in Broadhurst et al. 2005b and 
Medezinski et al. 2008. 

3. MEASUREMENTS OF LENSING EFFECTS 

These data permit high quaUty weak lensing measurements, 
following estabUshed methods required to deal with instru- 
mental and atmospheric effects following the formalism out- 
lined in Kaiser, Squires, & Broadhurst (1995), with modifi- 
cations described in Erben et al. (2001). The mean residual 
stellar ellipticity after PSF correction is less than, ~ 10"^, in 
all fields, averaged over 400-600 stars and consistent with the 
standard error on this measurement, ~ 10""*. Note also the 
lack of any systematic deviation from zero in the B-mode, gx, 
shown in Figure 2a. Full details of the methods are presented 
in Umetsu & Broadhurst (2008). 

Figure 1 shows the gravitational shear field by locally aver- 
aging the corrected distortions of colour-selected background 
galaxies of each cluster This is compared to the distribution 
of color-selected cluster sequence galaxies. In each case, one 
large syrmnetric cluster is visible around which the lensing 
distortion pattern is clearly tangential, with little significant 
substructure. The derived radial profiles of the lensing distor- 
tion are seen to be very similar in form (left panel of Figure 
2), with differing amplitudes reflecting a range of mass. 



The NFW profile fits well each cluster (Figure 2; Table 1), 
but surprisingly the derived concentrations lie well in excess 
of the standard ACDM model (Figure 3). For the derived 
masses of our clusters, the predicted mean concentration is 
Cvir ~ 5.06 ± 1.1 based on Duffy et al. (2007) using the im- 
proved WMAP5 (Komatsu et al. 2008). Our best-fitting val- 
ues instead range over 8 < Cvir < 15, with a mean value of 
((cvir) = 10.39 ± 0.91) representing a 4.8(t discrepancy. Our 
measurements also imply A370 is the most massive cluster 
now known, Myir = 2.93 x lO^^MQ/h^o, with a virial mass 
twice that of RXJ1347-11, the most X-ray luminous clus- 
ter known. Our estimate of the virial mass of RXJ1347-11 
is in very good agreement with independent X-ray and lens- 
ing analyses (Miranda et al. 2007, Kling et al. 2005). For 
A1689, a somewhat lower concentration parameter is derived 
by Limousin et al. (2007), Cvir = 9.6 ±2.0, from indepen- 
dent weak lensing measurements, which is consistent with our 
findings at the l.Sa level, whereas the observed Einstein ra- 
dius of approximately 52"is under-predicted by the NFW pa- 
rameters obtained by Limousin et al. (2007), 24" ± 1 1" (see 
Umetsu & Broadhurst 2008). 

Note, when model fitting an estimate of the background 
depth is required. The mean depth is sufficient for our pur- 
poses as the variation of the lens distance ratio, Dds/Dg, is 
slow for our sample because the clusters are at relatively low 
redshift compared to the redshift range of the background 
galaxies. The estimated mean depth of the combined red-nblue 
background galaxies is listed in column 5, which we obtain by 
applying our colour-magnitude selection to Subaru imaging 
of the HDF region (Capak et al. 2004), where photometric 
redshifts are reUable for red galaxies, with a mean redshift 
close to z = 1.2 depending on the details of the colour mag- 
nitude selection. For the blue galaxies we may rely on the 
zCOSMOS deep redshift survey (Lilly et al. 2008) where our 
mean depth is typically z — l.l. 

Lensing measures projected mass and so a statistical bias 
arises from the triaxiality of clusters in cases where the ma- 
jor axis lies along the line-of-sight. This leads to an ^ 18% 
increase in the mean value of lensing derived concentrations 
based on ACDM (Oguri et al. 2005; Hennawi et al. 2007), 
and an overall discrepancy in c™ of 4.0(7 with respect to the 
predictions. A larger bias is inferred for clusters selected by 
the presence of large arcs, ^ 34%, representing the most triax- 
ial cases (Hennawi et al. 2007). Our sample is defined by the 
quahty of available imaging and includes clusters observed 
for reasons other than lensing, hence it is unUkely that these 
clusters are all particularly triaxial with the long axes pointing 
to the observer Even so, applying the maximum estimated 
bias (^ 50%, Oguri & Blandford 2008) cannot account for 
our measurements (Figure 3). 

Multiply-lensed images are visible in all our clusters, from 
which the inner mass distribution may be determined (Broad- 
hurst et al. 2005a; Halkola et al. 2008; Hennawi et al. 2008) 
and an equivalent Einstein radius derived by averaging az- 
imuthally (Figure 2) which for most cases is close to the ob- 
served radius. For A1703 and A370, the tangential critical 
curves are significantly elongated, for r < 1' (Limousin et al. 
2008; Kneib et al. 1993) although in general the mass distri- 
bution is always less elliptical than the critical curve. X-ray 
observations of A370 show the inner region has some sub- 
structure (eg. Shu etal 2008), but beyond r > 1' the weak 
lensing pattern and the distribution of galaxies is symmetrical, 
as seen in Figure 1. For each cluster the NFW parameters de- 
rived with or without the Einstein radius constraint are found 
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Table 1. the Subaru distortion measurements combined with the 
Einstein-radius constraint. 
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Cvir 




A1689 
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A370 
RXJ1347-11 
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0.258 
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33 fe, = 2.8) 
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0.066 


i.J5'_Q 22 

1 30+* 24 

2 03+0.36 
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7 75+'-'2 
-0,92 

^"•^ -2.13 


4.94/9 
2.69/5 
5.54/8 
6.25/7 



to be closely consistent and in combination they are more pre- 
cisely determined - see Table 1. 

We also examine the magnification profile, via the 
surface number density of background galaxies (Broadhurst, 
Taylor, & Peacock 1995). At faint fluxes where the counts fol- 
low a power-law slope, s = d\ogN{< m)/dm, lensing modifies 
the true density above the magnitude hmit , A'o(w < wiim), by, 
N{< miim) = ^o(< 'Wiim)M('')^'^'~^. implying competition be- 
tween the magnified sky, which reduces the surface density, 
and an increase of galaxies magnified above the flux limit. 
Here we use only galaxies lying redward of the cluster se- 
quence because for these the intrinsic count-slope of faint red 
galaxies is relatively flat, s ~ 0.1, so a net count depletion re- 
sults, which we readily detect for each cluster, increasing to- 
wards the center (Figure 2). For each cluster we calculate the 
depletion of the counts expected for the the best fitting NFW 
profile derived from the corresponding distortion measure- 
ments above (Figure 2, right panel, with same colour code), 
finding clear consistency, which considerably strengthens our 
conclusions, and also establishes the utility of the background 
red galaxies for measuring magnification. The variation in 
the total number of background counts seen in Figure 2 sim- 
ply reflects the relative depth of the imaging data and also 
to some extent the redshift of the cluster, which Ues redder in 
colour space with increasing cluster redshift, thereby reducing 
the numbers of background galaxies which can be selected to 
lie redward of the sequence for our purposes. Note, absolute 
counts are not used here, the models are instead normalised to 
the observed density when comparing with the model to avoid 
reliance on control fields. 

4. DISCUSSION 

We have shown that the mass distribution of a sample of 
the most massive clusters foUow a very similar form in terms 
of the observed lensing based profiles. The distortion pro- 
files measured are amongst the most accurate constructed to 
date and the great depth of the Subaru imaging permits mag- 
nification profiles to be established with unprecedented detail 
using the background red galaxy counts. The magnification 
and distortion profiles are fitted weU by the general NFW pro- 
file for CDM dominated halos, which continuously steepens 
with radius, but the profiles derived are more concentrated 
than predicted by ACDM based simulations for which clus- 
ters of ~ IO'^Mq, form late (z < 0.5) with low concentrations, 
reflecting the low mean cosmological mass density today. 

Tension with the standard ACDM model is also indicated 



by the observed Einstein radii of very massive clusters, in- 
cluding those studied here, which he well beyond the pre- 
dicted range derived from massive halos generated in ACDM 
based simulations (Broadhurst & Barkana 2008). In addition, 
a significant discrepancy is also claimed for X-ray derived 
cluster concentrations, with respect to the ACDM based on 
WMAP5 parameters (Duffy et al. 2008). The combined X- 
ray and lensing based concentrations compiled by Comerford 
and Natarajan (2008) also exceeds somewhat the standard ex- 
pectation, despite including X-ray based concentrations that 
are found to be systematically lower than for lensing (Com- 
erford and Natarajan 2008). The mass-concentration relation 
of Mandelbaum etal (2008) is claimed to be more consistent 
with ACDM. Here allowance is made for the smoothing of 
the central lensing signal when stacking many SDSS clusters, 
relying on the brightest cluster galaxy to define the center of 
mass. Cluster members misclassified as background galaxies 
when estimating photometric redshifts will also act to reduce 
the derived concentrations, diluting the central lensing signal. 

Without resorting to radical proposals regarding the nature 
of DM, our findings imply the central region of clusters col- 
lapsed earlier than expected. Assuming the simple redshift 
relation Cvir(z) = co(l +z) (e.g., Wechsler et al. 2006), relat- 
ing central cluster densities to the cosmological mean density, 
then the formation of clusters with Cvir(l+z) = 10-15 cor- 
responds to z > 1, significantly earlier than in the standard 
ACDM, for which clusters form today with cq ~ 5.6. 

The presence of massive clusters at high redshift (z ^ 2), 
and the old ages of their member galaxies (Zirm et al. 2008; 
Blakeslee et al. 2003), may also imply clusters collapsed at 
relatively early times(Mathis, Diego, & Silk 2004), for which 
accelerated growth factors have been proposed, adopting a 
generahsed equation of state (Sadeh & Rephaeh 2008). Al- 
ternatively, since clusters correspond to rare density maxima, 
then any non-Gaussianity in the early fluctuation spectrum 
may advance cluster formation (Matthis, Diego, & Silk 2004; 
Sadeh, Rephaeli, & Silk 2007). Our results present a chal- 
lenge to some of these models, since earlier cluster formation 
must not also enhance significantly the abundance of clusters. 
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Fig. 1. — Maps of the surface number-density distribution of color-selected cluster member galaxies, with the gravitational shear of background galaxies 
overlayed; 10% ellipticity is indicated top right, and the resolution of the distoilion map is shown bottom right. In each case a single concentration of galaxies is 
visible, around which a coherent tangential pattern is centered, with little significant substructure. 
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Fig. 2. — (Lefty. Tangential distortion profiles of the 4 clusters based on the combined red and blue background samples and compared with the best-fitting 
NFW model. The Einstein-radius determined from multiply-lensed images is indicated, marking the point of maximum distortion, g+ = 1. Shown below is 
the rotated shear, gx, demonstrating no obvious non-tangential distortion. (Right): Profiles of background red galaxy counts, whose intrinsic slope is relatively 
shallow, so lens magnification reduces the observed numbers towards the center Over-plotted are normalised NFW models derived from the distortion profiles in 
the left panel, demonstrating consistency between these two independent lensing observables; note A370 is shifted upward 40% for clarity. The count uncertainty 
is obtained by sub-dividing each annulus into equal area cells, with a tail of > 2a cells excluded to remove inherent small scale clustering of the background. 
Small areas around each bright object (mn < 21) are excluded, including cluster members, interior to 3 times the isophotal radius, where the detection of faint 
galaxies is significantly compromised (Broadhurst et al. 2005b; Umetsu & Broadhurst 2008). 
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Fig. 3. — Comparison of observations with the ACDM model, based on A'-body simulations for the Cvir~^vir relation, derived at z = 0. The predictions of 
Duffy et al. (2008) (erg = 0.8, WMAP5) and Neto et al. (2007) (erg = 0.9, WMAPl) are shown as solid curves, with Icr uncertainty (from Neto et al. 2007) 
indicated by the hatched area. Also shown is the level of selection and orientation bias for projected masses based on ACDM (dashed curves, see text). The 
data points are all derived from lensing alone, and are multiplied by (l+z)^'^^ at the cluster redshift, for consistency with the evolution of fvir(^vir) derived from 
ACDM simulations by Duffy et al. 2008. Note, a significant correction for substructure is required for C10024-I-16 (Umetsu et al., in preparation). Clearly the 
data he well above the predicted relation even when possible sources of bias are considered. 



